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ON  THE  HITCHCOCK  DISTRIBUTION  PROBLEM 

Merrill  M,  Flood 

1 .  Introduction 

Frank  L.  Hitchcock  [l]  1  has  offered  a  mathematical 
formulation  of  the  problem  of  determining  the  most  economical 
manner  of  distribution  of  a  product  from  several  sources  of 
supply  to  numerous  localities  of  use.  He  also  suggests  a  com¬ 
putational  procedure  for  obtaining  a  solution  of  his  system  in 
any  particular  case.  L.  Kantorovitch  [2] ,  Trailing  C,  Koopmans 
[^] ,  George  B.  Dantzlg  [4b] ,  C.  3.  Tompkins  [5] ,  Julia  Robinson 
[b]  ,  Alex  Order*  [7]  and  others  [4]  have  also  discussed  the 
computational  aspects  of  this  problem. 

The  present  paper  is  concerned  only  with  the  mathematical 
justification  of  computational  procedure ,  and  is  limited  to  one 
specific  method  of  solution  of  general  validity.  No  attempt  is 
made  to  compare  the  various  methods  already  proposed,  either  as 
to  their  mathematical  similarity  or  as  to  their  relative  effici¬ 
ency  in  any  particular  case.1 
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2.  The  Problem 


Th®  problem  la  to  find  &  set  of  values  for  the  urn  variables 


X1J’ 

subject  to 

the  following  conditions: 

4B 

n 

(2.1) 

1  X11  m 
1*1 

cr  '  ri' 

(2.2) 

XU  *  °> 

(2.3) 

.Vi.idu 

1  >  J 

•  minimum. 

Th*  BUMb«r»  a,  n,  ?1 ,  e«,  sat  are  oivw  positive  integers* 

with  £c^  -  Zr^ .  Th®  Indices  1  and  J  ar®  understood  always  to 

range  over  these  sane  integers.  It  la  also  assumed,  for  convenience, 

that  m  >  n.  Any  set  of  values  x.  ,  that  satisfies  ail  these  condi- 

-*■  J 

tleas  is  called  a  solution  of  the  preplan. 

X*  will  se»etiaea  he  sere  convenient  to  use  sa  alternative 
statement  of  the  problem,  In  matrix  notation,  as  follows: 

(2.4)  M*y  £  b, 

(2.5)  7  >  0, 

(2.6)  a'y  -  minimum. 

It  is  easily  seen  that  the  two  formulations  are  equivalent  If  y,  a, 
b,  and  M’  are  defined  as  follows: 

7n(i-i)+J  *  *1J’ 
an(l-i )+J  ‘  dlJ’ 
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I  I  •••  I 

n  n  n 
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I 

H* 

— c 

n  n  n 

b  « 

J4  Ja  *•*  Jm 

r 

-~Ji  *  *  *  — Jm 

-r 

w her®  I  Is  the  Identity  matrix  of  order  n,  and  J.  is  the  mxn  matrix 
n  l 

with  all  elements  zero  except  for  the  i—  row  in  ¥hich  each  element 
is  unity.  Of  course,  y,  a,  c,  and  r  are  column  matrices  (or  vectors) 
with  components  an(l-j)+J’  Cj  and  ri ’  rt®P*ctlvoly* and  m 

prime  denotes  the  transpose  of  a  matrix  for  vector). 

3*  Fundamental  Theorems 

There  are  several  fundamental  theorems  concerning  systems  of 

linear  inequalities  that  are  useful  for  this  paper.  I  reproduce 
their  statements  here  in  a  form  due  to  A.  W.  Tucker.4  The  interested 
reader  can  find  proofs  of  these  theorems,  and  of  others  of  similar 
type,  In  a  paper  by  Gale,  Kuhn,  and  Tucker  [4c]. 

Fundamental  Problems;  (Here  lower  case  roman  letters  denote 
one  column  vectors,  while  capitals  denote  rectangular  matrices;  H#, 
a,  and  b  are  given  but  d  is  to  be  determined.) 

e 

Problem  1.  To  satisfy  tbs  eonstrelnts  *x  <  e,  x  £  0,  end  make 

b'x  -  d  for  d  maximal  in  the  sense  that  no  x  satisfying 
the  constraints  makes  b'x  >  d. 

Problem  II.  To  satisfy  the  constraints  M’y  £  b,  y  £  0,  and  make 

a'y  -  d  for  d  minimal  In  the  sense  that  no  y  satisfying 
the  constraints  makes  a'y  <  d. 


Problems  I  and  II  are  said  to  be  dual. 


Fundamental  Feasibility  Theorem: 

The  constraints  in  a  problem  are  feasible  (i.e.,  satisfied  by 
some  x  or  y)  If  and  only  if  the  dual  problem  in  homogeneous  form 
(i.e.,  with  b*0  or  a**0;  ns**  a  null  solution. 

Fundamental  Existence  Theorem: 

i.  x  and  y  are  solutions  of  Problems  I  and  IT  if  and  only  if 
they  satisfy  their  constraints  in  the  two  problems  and  uiatte  a’y  -  b 
Such  x  and  y  exist  if  the  constraints  in  both  problems  are  feasible 
ii.  A  problem  has  a  solution  if  and  only  if  its  constraints 
are  feasible  and  its  homogeneous  form  has  a  null  solution. 

Fundamental  Duality  Theorem: 

A  problem  haa  a  solution  (for  a  unique  d)  if  and  only  if  the 
dual  problem  has  a  solution  (for  ths  same  d). 

4.  The  Dual  and  Combined  Problems 
We  note  that  the  problem,  as  stated  in  relations  (2.4)-(2.6), 
is  a  fundamental  problem  of  form  II.  The  dual  problem  is; 

(4.1)  Hx  <  a, 

(4.2)  x  £  0, 

(4.3)  b'x  -  maximum. 

This  can  be  rewritten  in  a  more  convenient  form,  for  our  present 
purposes,  as  follows: 

(*.4)  Vj-Uj  <  dir 

(4.5)  ^jvj  ~  iriui  •  maximum, 

where  v.  -  x.~x  .  , ,  and  u.  -  -x_  . ♦x__ . 

J  j  n+j  1  *n+l  tn+m-n 


Theorem  1.  The  problem  has  a  solution. 
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Proof:  By  the  Fundamental  Existence  Theorem , there  la  a 
solution  If  and  only  If  the  constraints  are  feasible.  j.nd  >  •  0  is 

a  solution  of  the  problem  wtw  b  *  0.  He*  •  Sr^  #  *o 

X1J  “  **  a®®*****®*®*  ^ 

the  only  values  that  satisfy  the  constraints  are  x^j  *  0,  and  so 

the  theorem  is  proved. 

By  the  Fundamental  Duality  Theorem,  we  see 
Corollary  1A.  The  dual  problem  has  a  solution. 

Theorem  2.  The  numbers  x^,  and  ui,vj'  are  solutions  of  the 

problem  tad  the  dual,  respectively,  If  iM  «aly  If  tWbjr  eefcieTye 

<4-6>  ixjj  -  r,,  at.j-Oj. 

(..7)  >  0, 

(4.8)  Xlj(dlj+'J1-Vj)  ' 


Proof;  Since  (4.6)  and  (4.7)  are  simply  the  constraints  for 
the  problem  and  the  dual,  respectively,  it  remains  only  to  show 
that  (4.8)  Is  equivalent  to  the  condition  a’y-b'x  -  0.  Now 

A ' y— b ' x  -  f /j  +  flul 

'i!jXiJdij  "i  jXijVj  +ijIijUi 


Since  each  terra  in  this  sum  is  non-negati ve ,  a'y— b’x  •  0  if  and 
only  If  ^(d^V-V  -  °' 


r-2. 


We  refer  to  the  oroblem  of  finding  values  for  x.  .,'u,,  and  v, 

1  J  1  J 

tnat  satisfy  (4 . d)-(4 . c)  a3  the  "combined  problem",  and  note  that 
toe  combined  problem  always  has  a  colit'on. 


5 ♦  Linear  Graphs 

It  win  be  convenient,  for  some  purposes,  to  associate  linear 
graphs  [u]  wltn  certain  subsets  of  the  elements  of  a  matrix 
S  ■  ||s  if.  If  I  is  a  given  subset  of  the  elements  of  S,  we  define 
the  I— graph  L  of  S  as  follows:  the  vertices  of  L  are  all  the  points 
(h,k)  in  the  Cartesian  plane  for  wnlcn  s^  £  I;  the  arcs  of  L  are 
aj.1  line  segments  Joining  pairs  of  neigh  coring  vertices  with  either 
eq *a±  abscissas  or  eq  iai  ordinates,  where  two  vertices  with  equal 
abscissas  (ordinates)  are  neighboring  if  they  are  not  separated  by 
another  vertex  with  the  same  absciss  i  (ordinate).  For  the  moment, 
denote  the  vertices  of  L  by  symbols  a,t,c,***,f  and  the  arcs  by 
symbols  s <ch  as  ab ,  be,***,  cf  (no  distinction  is  made  between  the 
arcs  ab  ana  -a)  a  Then  a  chair,  is  a  set  of  one  or  more  d'stinct 
arco  that  can  be  arranged  as  ab,  be,***,  de,  ef,  where  vertices 
denote!  \y  different  symbols  are  distinct.  A  eye  :e  is  a  set  or 
•distinct  ircn  (at  least  four  are  necessary)  thr.t  can  be  ordered  as 
at,  be,"*,  ef,  ft,  the  vertices  being  distinct  «b  in  the  case  of 
a  cnaln.  A  trtph  la  connected  if  each  pa^r  of  vertices  Is  joined 
by  a  chain.  *  ; i s t  is  a  graph  containing  no  cycles,  and  a  tree 
le  a  connected  forest. 

If  L  contains  v  ertices,  a  arcs,  and  p  connected  pieces, 

the  r.  ..roer^  *  a-  /  +  •  ’  s  Known  as  the  eye ,  o.motlc  nurd  er  (or  first 


P-213 

-7- 


Betti  number)  of  L.  It  follows  from  a  well  known  theorem  [8] 
concerning  linear  graphs  In  general  that:  (l)  L  1b  a  forest  If  and 
only  if  ><-  0,  and  ( 1 1 )  L  contains  just  one  cycle  if  and  only  if 
M  *  1. 

Note  that  L  contains  a  cycle  If  and  only  If  there  is  a  subset 
of  I  that  can  be  arranged  as  a  sequence 

8hikt'  “hjka*  sh2k8'  8h*k3'  *  *  *  ^h^'  shaki 

where  the  h's  and  k*s  are  distinct  among  themselves;  and  L  contains 

a  single  cycle  If  and  only  if  I  contains  Just  one  *.uoset  that  can 

be  arranged  In  tne  displayed  form.  We  call  such  a  subset  of  I  an 

I — c 1 rc ul t  on  S,  and  denote  It  by  [s^] .  For  a  particular  arrangement 

of  [S J ,  we  shall  also  refer  to  the  ter~s  s.  .  as  odd-terms,  the 
<T  h^  - 

others  as  even-terms. 

In  case  I  consists  of  all  shlc  >  0,  as  It  frequently  will,  we 
speak  of  the  positive  graph  of  S,  positive  circuits  on  S,  and 
abbreviate  s  ich  statements  as  "the  positive  graph  of  S  is  a  forest" 
to  "S  1 s  a  forest". 

6.  Tne  Method  of  Solution 

In  the  met.uod  of  solution  to  be  developed  for  the  problem, 
we  start  with  a  special  set  of  values  X°  -  that  satisfy  the 

constraints  (4.6).  '.Ve  then  test  to  determine  w*' ether  or  not  there 
exist  and  v^  satisfying  the  re.ations  (**.7}  and  (4.8)  for  the 
plven  X°.  If  sc,  then  X°  Is  a  solution,  otherwise  not.  The  method 
next  ylelus  n  new  trla.  matrix  X1  -  |xj,|,  if  X°  Is  not  a  solution, 

*  t.1 

such  that  .£,(x? .-x‘  )d.  .  >  ..  After  a  finite  number  of  steps  this 

*  J  *  v  *  J 
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procesi  necessarily  must  terminate  and  it  leads  to  an  exact  integral 
solution  of  the  problem. 

The  first  trial  matrix  X°  is  a  forest  of  t  trees,  and  has 

m+n-t  non-zero  elements.  According  as  t  ■  1  or  t  >  1,  two  essentially 

6 

different  cases  may  be  met  at  each  stage  of  th*  solution  process. 

At  each  stage  when  X  *  (xij|  i8  a  tree,  the  equations  (4.8) 

have  a  general  solution  for  and  v^  with  one  fre*.  parameter,  say 

Uj.  However,  the  quantities  d^  i+ui“vj  are  uniquely  determined  in 

this  case,  so  It  is  sufficient  to  calculate  them  and  note  whether 

or  not  they  are  all  non-negative  in  order  to  decide  whether  or  not 

X  is  a  sol  ition.  If  some  d,  .  +u.  -v.  <  0,  then  there  is  a  uniq.e 

mji  h  Jt 

I-clrcult  fx  1  on  X,  wnere  I  consists  of  x.  «  and  all  positive  x. 

‘  B  lljl  1J 

that  may  be  arranged  with  x.  .  as  the  seoond  term,  say.  Let 

1  J  I 

g  aenote  the  smallest  odd-tenn  of  [X.] .  Then  the  new  tria1  matrix 
X*  is  obtained  from  X  by  adding  g  to  the  even  terms  of  [Xfl] , 
subtracting  g  from  the  odd-terms,  and  leaving  the  other  elements 
of  X  unchanged. 

At  each  stage  when  X  Is  a  forest  of  t>±  trees,  the  eqiations 
(4.8)  have  a  general  sol  ition  fbr  and  v.  with  t  Independent 
parameters,  and  the  quantities  d^ ,+  ‘i“vj  invo*ve  t-i  Independent 
parameters.  The  rows  and  columns  of  the  matrix  X  are  rearranged  so 
that  it  can  be  represented  as  a  sq  iare  matrix  of  order  t  whose  t2 


eLeaents  are  a  »bma trices  XaK  such  that  Xt,  ■  0  If  a  /  t  and  XOQ 

&C  au  aa 

la  a  tree  with  m  +n  -1  non-zero  elements  and  Is  of  order  m  xn . 

a  a  a  a 

We  can  select  U| ,  u,  . .  to  be  the  t  parameters. 

If  we  assign  these  the  value  zero  and  denote  this  particular 
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aoiution  of  (4.6)  by  i,  and  v,  then  we  may  define  numbers 
Pj.  j  *  We  Part^tion  the  matrix  ?  •  |j ^ }J  into  submatrices 

corresponding  to  the  X&b  and  denote  ths.r.  Tab.  Let  be  the 
smallest  element  in  F  and  define  the  square  matrix  F  of  order  t 
by  P  •  i|pabH.  To  designate  the  position  of  paf  in  the  matrix 
?  *  || PA  ,||,  we  may  write  pab  alternatively  as  paab\  the  subscripts 
referring  to  the  s  .bmatrix  and  the  superscripts  to  the  rows  and 
columns  In  the  submatrlx. 

1 

The  test  ae  to  whether  or  not  X  is  a  aoiution  consists  of 


p  +P 

^aja a  va,a3 


♦F.  for 

ahai 


forming  all  sums  p 

a  1 a  a  *  * • a b 

n  -  2,3,“*,t  ,  wnere  (ala***«ar)  Is  any  permutation  of  h  different 
positive  integers,  none  greater  than  t;  X  is  a  solution  if  and  only 
if  all  such  s  -ms  are  non— negatl ve. 

If  any  <  0,  then  there  is  a  unique  I— clrc  it  [x  1 

on  X,  where  I  consists  of  ail  positive  x.  ,  together  with  all  x,  , 

U  v  ^  <j 

ia  Ja 

_  ®  lr  ®  j  ±  | 

thet  correspond  to  the  terms  p_  of  p  .  _  , 

!a  V 

which  can  bs  arranged  to  involve  a.i  x  *  *® 

*k  ?k+i 

even-terms.  If  g  is  be  smallest  od  i  terra  In  [l_] ,  then  (as  in 
the  norv-degenerate  case)  the  new  trla.  matrix  X*  is  obtained  by 

adding  g  to  the  even-terms  of  [X.l ,  s  .btracting  g  from  the  odd- 

O 

ter me,  and  leaving  the  other  elements  of  X  ^changed. 


7 .  Tr.e  Initial  Trla.  Sol  .tlon 
An  X  tnat  satisfies  (4.o)  wi*.  be  cai.ed  a  trl_a_i  so  i  .fen. 
It  wo  ;ud  be  a*±  rigr.t  to  take  tr.e  positive  vai.es  r,c,/Ir,  for 
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the  Initial  trial  solution  X°  •  }} j  II •  An  alternative  is  to  construct 
an  initial  trial  solution  that  is  a  forest.  It  is  always  possible 
to  do  this  in  integral  values.  The  following  theorem  certifies  the 
existence  of  such  an  integral  trie  1  solution.  The  method  of  proof 
shows  hew  to  construct  one. 

gfieorem  j.  There  Is  a  matrix  X°  -  Ijx^il  with  integral  elements 
that  satisflos  (4.6)  and  is  a  forest. 

Proof t  The  theorem  Is  trivial  for  a-1.  Assume  the  theorem  Is 

true  for  m  and  consider  the  case  m+1. 

Let  the  notation  be  chosen  so  that  ri  >  r*  >  • • •  >  r  s  o 

-  -  -  m+i  '  ’ 

and  Ci  >  ca  >  *  *  *  >  cn  >  0.  If  n  <  m+1,  then  cx  >  r  +|.  If  n  *  m+1 
then  ci  >  rm+l  unless  c1  -  r^  -  A (for  all  i  and  ft  in  this  latter 
case  x°  »  A  satisfies  the  cond’t^ns  of  tp*  theorem.  Hence,  by  the 

induction  hypothesis,  there  is  s  set  of  non-negative  integers  x* 

•  •  ij 

(1  “  !>•••,*)  such  that  Jx*  -  cj^ijrm^i  "  *%  ,  and  X*-  |xjj| 

is  a  forest.  Then  X°  defined  by  xjj  •  x^,  x£+j  ^  •  c*Jrm+i»  satisfies 

(4.6) .  Now  since  the  (m+1 ~ ' row ,  with  only  one  positive  element, 
clearly  cannot  contribute  tt^as  to  a  positive  circuit,  X°  is  also 

a  forest ;  and  the  theorem  is  provcij 

TO  apply  this  method,  in  the  construction  of  a  trial  solution, 

search  for  the  smallest  r.  and  the  largest  c,  ,  and  then  set 

1 1  Ji 

x?  .  -  r  .  In  effect,  this  deletes  the  ix—  row,  after  c  .  is 

1  1 J 1  1  ,]j 

replaced  by  c 1  -r.  ,  and  the  process  is  repeated  (interchanging 
rows  and  columns  at,  necessary’)  until  all  have  been  determined. 

For  a  itomatic  machine  calculation  the  procedure  is  easily  made 
unique,  for  any  one  starting  order  of  rows  and  columns,  by 


'  1-2 12 

-11- 

specifying  that  the  search  Is  first  on  row-totals  when  the  number 

of  rows  1 s  the  same  as  the  number  of  columns  at  any  stage,  and  that 

the  row-total  or  column-total  with  the  smallest  Index  is  chosen 

whenever  at  any  stage  there  are  several  equal  values  to  choose  from. 

This  initial  trial  solution  will  be  called  "preferred"  for 

7 

identification. 

Theorem  4.  A  trial  solution  that  is  a  forest  of  t  trees  has 
m+n-t  non-zero  elements. 

Proof :  Observe  first  that  if  the  trial  solution  X  is  a  forest 
of  t  trees,  the  rows  and  col  imns  of  X  can  be  rearranged  so  that  X 
has  the  form 

Xu  0  •  •  •  0 

0  X*a*  *  *  0 


0  0  •••  Xtt  II 

where  each  Xaa  is  a  tree.  Consequently,  the  theorem  amounts  to 
proving  that  an  mxn  matrix  with  no  zero  rows  or  columns ,  which  is  a 
tree,  nas  m+n— 1  positive  elements.  If  m+n  ■  2,  this  is  obvious,  so 
assume  the  statement  to  be  true  for  all  matrices  for  which  m+n  «  k  and 


consider  one  for  which  m+n  *  k+1.  Since  m  >  n,  clearly  some  row 
has  only  one  positive  element,  as  otherwise  there  would  be  a 
positive  circuit.  Delete  this  row  and  apply  the  induction  hypothesis. 
In  actual  cases  when  m  and  n  are  relatively  small,  or  when 


there  is  other  reason  to  believe  that  an  initial  trial  solution 
better  than  the  preferred  one  can  be  found  by  trial  and  error,  it 
may  be  better  to  constrict  the  initial  trial  solution  In  some  other 
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way  than  the  one  given  in  the  proof  of  Theorem  3,  in  order  to 
reduce  the  number  of  steps  required  in  the  iterative  process. 

The  methods  developed  in  this  paper  apply  directly  for  any 
trial  solution  that  is  a  forest,  and  are  readily  extended  for  other 
cases.  It  is  easy  to  see  that  there  must  be  at  least  one  solution 
which  is  a  forest. 


8.  Non-degenerate  Case 

We  consider  now  the  case  of  a  trial  solution  X  which  is  a 

tree.  Let  the  positive  elements  of  X  be  x,  ,  ,  a  *  1,»* • ,m+n-l. 

1a',a 

We  shall  need  the  following  theorem. 

Theorem  5.  If  X  Is  a  trial  tree,  the  set  of  equations 

(8*1)  diJ+ui~vJ  "  0  for  (i,J)  *  (1a,4^ 

has  the  general  solution 

ai  "  ui+z»  vj  "  vj+z' 

where  Is  a  particular  solution  and  z  is  arbitrary. 

Proof:  The  theorem  is  apparent  for  m  ■  1,  and  we  proceed  by 
induction..  Suppose  the  theorem  is  true  for  all  trial  trees  of  m 
rows,  and  let  X  be  an  m+lxn  trial  tree.  Obviously,  there  must  be 
at  least  one  row  of  X  that  has  exactly  one  non-zero  element;  we 
may  suppose  It  to  be  Xj^  without  loss  of  generality — also  that 
i  ■  m+1  and  J  .  *  n.  Since  X  Is  a  trial  tree,  the  matrix 

obtained  from  X  by  deleting  the  last  row  (or,  if  m+l  -  n,  its 
transpose)  is  also.  The  induction  hypothesis  implies  that  the 
general  solution  of  (8.1),  with  the  final  equation  omitted,  Is  of 
the  form  u^  -  uj+z,  v^  *  v*+z.  We  note  next  that  this  final  equation 
becomes  utn+  -  (vn-dm+in)  +  z  -  u^,  +  z.  The  theorem  follows  easily. 

It  will  oe  convenient  to  call  the  particular  solution  u^,  7^ 
of  (8.1)  obtained  by  setting  ut  -  0  the  preferred  trial  solution 
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of  the  dual  problem  corresponding  to  the  trial  tree  X.  As  an 
obvious  consequence  of  Theorem  5>  we  state 

Corollary  5A.  If  X  is  a  trial  tree,  then  it  is  a  solution  of 
the  problem  if  and  only  if  the  corresponding  preferred  trial 
solution  (u1#Vj)  of  the  dual  problem  satisfies  £  0 

(Tor  all  I  and  J).  J 

All  that  is  needed  now  in  order  to  establish  the  method  for 


the  non-degenerate  case  is  to  show  how  to  construct  a  new  trial 
matrix  X*,  if  X  is  not  a  solution,  such  that  2  (x.  .-x*  )d,  .  >  1. 


1J  i J  i J  - 

dk pvn  < 0  eor 


In  this  case,  it  follows  by  Corollary  SA  that 
at  least  one  pair  (k,/)  and,  of  course,  x^  ■  0. 

Theorem  6.  If  the  trial  solution  X  is  a  tree,  and  ■  0, 
then  there  is  a  unique  I-circuit  on  X,  where  I  consists  of  all 
positive  x^  together  with  x^. 

Proof:  It  suffices  to  show  that  the  I— graph  of  X  has  cyclomat I c 
number  1.  By  assumption,  the  positive  graph  of  X  has  cyclomatic 
n;inber  zero,  and  since  X  must  have  positive  elements  and  x^^ 
for  some  a  and  b,  the  I-graph  of  X  has  two  more  arcs,  one  more  vertex, 
and  the  same  number  (one)  of  connected  pieces.  Hence  1,  and 


the  proof  is  complete. 

Now  arrange  this  unique  I— circuit  [Xfl]  with  x,^  as  the  second 
term  and  let  g  be  the  minimum  of  the  odd-tenns  of  [Xg]  in  this 


arrangement.  If  we  subtract  g  from  the  odd-teme,  add  g  to  the 
even  terms,  and  leave  the  remaining  elements  of  X  unchanged,  we 
get  a  matrix  X  that  satisfies  (4,  j)  and  is  a  forest  (since 
was  unique). 
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Theorem  7.  S  (xJ  ,-x*  ,)d,  ,  >  1. 
^  j  ij  ij  l j 


Proof 


!  ut  -  D4,Jt.  *1,0,*  *i,o,'  *i,o,*—*  xV.'\j,] 


where  x.  .  -  x,  ,.  Then 

iiJt  kt 


1^J(*10"xio)dlJ  ■  s(dl,Ji”di,J,+<!i,0.‘<3i.Js+""f'1isJB~disJ,) 

■  -*<dl,0.+«l,-w0.)  *  x- 


The  theorem  follows. 

If  X#  Is  a  tree,  then  the  whole  process  is  repeated  until  at 

some  stage  a  trial  matrix  Is  obtained  that  either: (3)  is  a  solution, 
or(ii)  Is  not  a  solution  and  Is  a  forest  of  t  >  1  trees.  We  discuss 
(11)  next. 


9«  The  Degenerate  Case 

Let  X  be  a  trial  matrix  which  is  a  forest  of  t  >  1  trees.  As 
we  have  seen,  we  may  suppose  that  the  rows  and  columns  of  X  are 
ordered  so  that 

Xu  0  •  •  •  •  •  0 

o  Xaa  •••  o 


where  each  submatrix  Xaa  of  order  maxna  is  a  tree.  We  can  apply 
the  methods  of  the  non-degenerate  case  to  the  subproblems  corresponding 
to  the  submatrices  X  a  and  either  obtain  a  solution  to  ea  u  subproblen 
or  further  decompose  the  matrix  X,  so  we  may  also  assune  that  each 
Xaa  Is  a  solution  to  its  s ubproblem. 


By  Corollary  5A,  we  know  that 


(9.1)  4^  +  i*  -  7t*  >  °> 


l,-".ma;  Jg  -  !,-• 


where  u  .  v  Is  the  preferred  trial  solution  of  the  dual  subproblem 

a  fl 

corresponding  to  the  solution  Xa_,  and  that 

afil 


(9-2)  9aaaa  +  >  -  vaa  -  0  If  xaaaa  >  0. 


1. 


We  recall  also  that  the  moat  general  values  for  t?aa  and  vft*  are  given 
by 


'a  Ja  4  , 

u*  -  a  4  “a ' 


^a  -*^a 

y  »  V  -f  Z 

a  a  a 


where  the  2ft  are  arbitrary  parameters. 

It  follows  from  Theorem  2  that  X  Is  a  solution  If  and  only  If 
there  are  values  of  z  that  satisfy  Inequalities  corresponding  to 

a 

(4.7),  or  in  oar  present  notation: 


(9.3)  daabD  +  uaa  -  vbD  >  0  for  all  a,  b,  i#,  and  Jb. 

But  (9*3)  has  a  solution  for  z  If  and  only  If  the  following 

cl 

inequalities  have  a  solution  for  z  : 


(9.1*) 

where 


Pab+2a-zb  >  °- 


1  2 
-  a^b 

pa  b 


Vb 

AD  A 


p„b  b  j 


) 


4 


We  have  proved 

Lemma  A.  X  Is  a  solution  If  and  only  If  there  are  real 
n  imbers  z&  s  ich  that  Pat+za~zb  £  0,  a,b  -  l,**«,t. 


i 


-lo¬ 
in  order  to  establish  a  criterion  for  the  solvability  of  (9»^)» 
we  consider  a  special  case  of  the  original  problem,  defined  as 
follows:  dal3  »  pafe,  rfl  -  c^  •  1,  a,b  ■  l,***,t.  We  call  this  the 
special  problem,  the  corresponding  dual  the  special  dual ,  and  now 
consider  the  special  combined  problem: 

fab  *  fab  ‘  yab  £  °- 

Pab^a^b  *  °- 
yab^ab^a^b)  '  °* 

If  we  set  yab  •  then  for  this  trial  solution  the  conditions 

reduce  to: 

Pab^a^b  >  0  for  a  *  b> 

“  0> 

Since  paa  -  0,  it  follows  that  -  w&,  and  so  these  conditions 
are  equivalent  to  (9«^)«  Hence,  by  Theorem  2,  (9.^)  has  a  solution 
if  and  only  if||6aJ|is  a  solution  of  the  special  problem.  Using 
Lemma  A,  we  now  have 

Lemma  B.  X  is  a  sol.t'on  if  the  original  problem  if  and  only 

if  the  identity  matrix  is  a  solution  of  the  special  problem. 

Theorem  8.  X  is  a  solution  of  the  problem  if  and  only  if 
paia#...ah  £  °'  h  *  2,3,*««,t,  where  (at ,at,* • • ,ah)  is  any  permutation 
of  h  different  positive  Integers,  none  greater  than  t,  and 

pataa,#*ak  ^SiSt^*#***  *pa^ai* 

* 

Proof i  By  Lews  B,  it  suffices  to  short  that  the  ooadltlon  of 

the  theorem  is  equivalent  to  the  statement  that  l|6abll  is  a  solution 
of  the  special  problem. 

First  of  all.  It  is  easy  to  see  that  at  least  one  solution 
V  *  °**  the  special  problem  Is  a  forest,  and  hence  has  leas 

than  2t  non-zero  elements.  That  the  elements  of  Y  are  all  either 


t-'±} 

-17- 


zero  or  unity  can  be  seen  by  Induction  as  follows.  The  basis  of  the 

induction  is  obvious,  and  we  consider  the  case  t+1,  assuming  the 

statement  for  t.  There  must  be  at  least  one  element  of  Y  that  is 

unity,  as  otherwise  Y  would  have  at  least  2t  non-zero  elements.  We 

may  suppose  that  this  element  is  But  then  the  induction 

hypothesis  implies  that  each  element  y&b,  a,b  •  l,***,t,  is  zero  or 

one.  It  follows  that  there  are  exactly  t  elements  of  Y  that  are 

unity,  whence  L  y^.p,.  can  bo  written  as  p  ,  +  p  ,  +. ..+p„  v 
a,’o  ab  ab  a»bi  Ka»b*  atbt 

where  (aia2»**a|.)  and  (bib** ••b^)  are  pennutations  of  the  first  t 
integers.  Then  |j5abfl  1®  a  solution  of  the  special  problem  if  and 
only  if  always 

palbl+pa,b,+"-+patbt  >  Pi *+P*«+-  •  -+Ptt  *  0. 

Tne  proof  is  completed  by  noting  that  this  sum  can  be  written  as 
o  +P.  q  (aia*«**ah)  a®  described  in  the  theorem. 

d  J  eg  Cl  3  44 

We  now  need  to  show  how  to  construct  an  improved  trial  solution 
X*  in  the  event  that  X  is  not  a  sol  -tlon.  In  this  case,  *e  know 
from  Theorem  o  that  there  is  a  sum 


1°  ,0  jO  .0  jO  to 

_  ajht  _  a#Jaa  Va» 

n  +  n  +  •  •  •  •¥  D  <  0. 

*3i  a2  *aa  a2  ai 

iaJa 

►  .  _  _  j  _  i.  _  a  _  ^  _ j  a,  a _  i  ^  .  _  ** 


Let  I  consist  of  all  positive  elements  xa 
1° 

x„  *  of  X.  Then  we  assert; 

ak  ak+i 


together  with  all 


Theorem  9.  There  is  a  unique  I-clrcuit  on  X  that  can  be 

i°  j° 

Vak+» 

arranged  to  involve  sll  tne  x  as  even-terms. 

ak  ak+i 


Proof:  The  positive  graph  of  X  has  m+n-t  vertices,  m+n-2t  area, 

al# 

and  t  connected  pieces.  Also,  for  each  x  there  are  non-zero 

■it'll  |i|®  ^  k+: 

Vak  .  ak+.  \+i 


elements  x 


ak  ak 


»  x 


ak+i  ak+» 


Hence  in  passing  from  the  positive 


graph  to  the  I-gfaph ,  h  vertices  and  2h  area  are  added,  and  the  number 

of  connected  pieoea  la  decreased  from  t  to  t-h+1.  Thua  the  cyclo- 

aatlc  number  of  the  I-graph  is 

«  (2h-Hn+n— 2t)  -  (h-Hn+n-t)  +  (t-h+l)  »  1, 

so  there  la  a  unique  I-circuit  [X_]  on  X.  Since  the  graph  obtained 

i°  j° 

VVn  r  -1 

by  omitting  from  I  any  x  _  clearly  has  no  cycle,  (X  J  contains 

ak  0 

all  of  these. 

Evidently  [xj  can  be  arranged,  for  example,  as 


o  ,o 


i:  j.  i:  jr  i.  j: 

a j  at  at  a 2  a* 

:a»  at  *  xat  a*  9  xaM  a,  * 


1? 

a*  a# 


i!  J l 

•a 


•  •  •  ,  x 


a*  a*  *  a8  a3 


1°  J° 

ak  ak+i 

so  taat  all  x.  _  appear  as  even-terms. 
ak  ak+i 


As  in  the  non-degenerate  cc*se,  let  g  be  the  smallest  odd-term 
in  [xj  (hence  g  >  0)  and  define  a  new  trial  matrix  X*  by  replacing 
the  elements  of  X  that  appear  in  [xj  by  new  ones  Increased  by  g 
for  even-terms  and  decreased  by  g  for  odd  terns;  the  other  elements 
of  X  are  left  unchanged.  Again  X*  satisfies  the  conditions  for  a 
trial  matrix.  To  complete  the  discussion  of  the  degenerate  case,  it 
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Proof:  Since  X  and  X*  differ  only  on 

[Xal  '  xi.J»’"'’  xv.-’ 


xl 

1 


f 


then 


The  proof  ie  completed  by  noting  that  d^j  ■  4,1(1 

hj  m  0  lf  X1J  >  °-  80  that  '  -«(p.,a,...a 


). 


me 
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F-jotnotes 


1.  Numbers  In  brackets  refer  to  bibliography  at  the  end  of  the  paper. 

2.  My  interest  in  the  problem  was  aro used  by  papers  on  transportation 
theory  presented  by  Koopmans  P>aj  and  Dantzig  [4bJ  at  a  conference  on 
linear  programming  in  Chicago  during  June  1949»  under  the  auspices  of 
the  Cowles  Commission  for  Research  in  Economics  of  the  University  of 
Chicago.  Several  other  papers  presented  at  this  conference  are  of 
closely  related  interest.  Professor  Koopmans,  in  his  Introduction 

to  the  Conference  Proceedings  [4],  also  discusses  the  background  and 
interrelationship  of  the  conferei.ee  papers — lncl  udlng  the  bearing  of 
some  of  these  on  the  Hi ter. cock  distribution  problem.  The  result!  of 
the  present  paper  have  been  presented  in  three  seminar  lectures: 
once  in  December  1949  at  The  RAND  Corporation  in- Santa  Monica,  once 
in  July  1950  at  the  Institute  for  Nanerical  Analysis  of  the  National 
Bureau  of  Standards  in  los  Angeles,  and  once  in  June  1951  at  the 
National  Bureau  of  Standards  in  Washington,  D.  C. 

I  am  especially  indebted  to  Dr.  D.  R.  Pulkerson,  who  has  given 
real  assistance  in  simplifying  notation  and  proofs  of  theorems,  for 
a  careful  reading  of  the  manuscript. 

3«  Th*re  is  no  loss  of  generality  in  assuming  that  the  d. ,  are 
positive  integers,  rather  than  rational  numbers,  since  theJprotlem  is 
essentially  unchanged  if  d. .  is  replaced  by  ad..+b  where  a  and  t  are 
any  positive  rational  msnbirs.  I  have  not  exa/tifned  the  case  in 
which  some  of  tne  quantities  r. ,  c.,  and  d. .  are  Irrational.  The 
only  effect  of  irrationality  on  the  res..ltBJof  the  present  paper  is 
a  possible  lack  of  convergence  of  the  Iterative  process  of  sol  tion. 
These  considerations  are  not  of  Importance  in  the  isual  applications. 


4.  Unpublished  note,  December  1949. 

5.  We  omit  the  condition  (4.2),  that  x  >  0,  alnce  this  imposes  no 
limitation  on  u,  and  v<. 

6.  These  are  the  non-degenerate  and  degenerate  cases  in  the  work  of 
D*ntzig  [4b] .  I  shall  »ae  these  terms  also.  The  method  of  solution 
developed  by  Dantzig  [4b]  for  the  non-degenerate  case  Is  essentially 
the  same  as  the  one  in  tnis  paper,  although  the  derivations  of  the 
results  are  quite  different.  Orden  [?]  has  subsequently  given  en 
elegant  metnod  for  reducing  the  degenerate  case  to  the  non-degenerate 
one,  as  an  extension  of  the  a-method  proposed  by  Danzig  L4bJ .  I 
believe  that  the  treatment  of  the  degenerate  case  provide*.  the  only- 
res  ults  in  the  present  paper  that  are  new,  or  at  least  fresh  for  the 
Hitchcock  problem,  and  also  of  some  mathematical  interest.  It  also 
seems  likely  tnat  the  metnod  given  here  will  often  be  more  efficient 
computationally,  in  the  degenerate  case,  than  the  Dantzig— Orden 

€ -method. 


7.  Sometimes,  as  in  t  is  Instance,  I  indicate  how  to  make  a  un.q.e 
choice  among  possible  alternatives  ot  each  computational  3v,ep  but 
usually  I  do  not.  It  is  necessary  to  do  t..is  In  order  completely  to 
routinize  the  comp -ting  steps,  of  co  rse,  bit  the  matter  presents  r.c 
difficulty  and  I  on  t  it  here. 
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